Interacting fermions in quasi-one-dimensional harmonic traps 
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Quasi-one-dimensional (quasi-ld) two-component Fermi gases with effectively attractive and re- 
pulsive interactions are characterized for arbitrary interaction strength. The ground-state proper- 
ties of the gas confined in highly elongated harmonic traps are determined within the local density 
approximation. For strong attractive effective interactions the existence of a molecular Tonks- 
Girardeau gas is predicted. The frequency of the lowest breathing mode is calculated as a function 
of the coupling strength for both attractive and repulsive interactions. 

PACS numbers: 



The study of cold quasi-ld atomic quantum gases 
presents a very active area of research. So far, most of 
the experimental and theoretical IE ^ IE ly| investi- 
gations have been devoted to quasi-ld Bose gases and, in 
particular, to the strongly-interacting Tonks-Girardeau 
gas, which can be mapped to a gas of non-interacting 
fermions 0j0, 111- Quasi-ld atomic Fermi gases have not 
been realized experimentally yet. 

The role of interactions in quasi-ld atomic Fermi gases 
has been studied mainly in connection with Luttinger 
liquid theory 0, 0|. Recati et al. investigate the 
properties of a two-component Fermi gas with repul- 
sive interspecies interactions confined in highly-elongated 
harmonic traps. In the limit of weak and strong cou- 
pling these authors relate the parameters of the Lut- 
tinger Hamiltonian, which describe the low-energy prop- 
erties of the gas, to the microscopic parameters of the 
system. The Luttinger model is used to analyze the man- 
ifestations of the uncoupled dynamics of spin and density 
waves (spin-charge separation). 

This Letter investigates the properties of inhomoge- 
neous quasi-ld two-component Fermi gases under har- 
monic confinement with attractive and repulsive inter- 
species interactions. The present study is based on the 
exact equation of state of a homo gene ous Id system of 
fermions with zero-range attractive |lll Il2l | and repulsive 
(T3l | interactions treated within the local density approx- 
imation (LDA). We calculate the energy per particle, the 
size of the cloud, and the frequency of the lowest com- 
pressional mode as a function of the effective Id coupling 
constant, including infinitely strong attractive and repul- 
sive interactions. Moreover, for attractive interactions 
we discuss the cross-over from the weak- to the strong- 
coupling regime and point out the possibility of forming 
a mechanically stable molecular Tonks-Girardeau gas. 

Quasi-ld two-component Fermi gases with effectively 
attractive and repulsive Id interspecies interactions can 
be realized in highly-elongated traps. The behavior of 
atomic gases tightly-confined in two directions can, if 
the confinement is chosen properly, be characterized to a 
very good approximation by an effective Id coupling con- 
stant, gld^ which encapsulates the atom-atom interaction 



strength. This coupling constant can be tuned to essen- 
tially any value, including zero and ±cx3, by varying the 
3d s-wave scattering length 03^ through application of 
an external magnetic field in the proximity of a Feshbach 
resonance. 

In homogeneous Id Fermi gases with attractive inter- 
actions, sound waves propagate with a well defined ve- 
locity, while spin waves exhibit a gap jl2j| . Furthermore, 
in the strong-coupling regime, the ground state is com- 
prised of bosonic molecules, consisting of two fermions 
with different spin whose spatial size is much smaller 
than the average intermolecular distance Conse- 
quently, BCS-type equations have been discussed for ef- 
fectively attractive Id interactions The quasi-ld 
molecular Bose gas discussed here has similarities with 
the formation of a molecular Bose-Einstein condensate 
(BEG) from a 3d Fermi sea close to a magnetic atom- 
atom Feshbach resonance llal , which may allow study of 
the BCS-BEC cross-over [l^. In contrast to the 3d case, 
the cross-over from the weak-coupling regime of large, 
overlapping Gooper pairs to the strong-coupling regime 
of small, tightly bound Bose pairs can be investigated in 
Id within an exactly solvable model. 

Consider a two-component atomic Fermi gas confined 
in a highly-elongated trap. The fermionic atoms are as- 
sumed to belong to the same atomic species, that is, to 
have the same mass m, but to be trapped in different hy- 
perfine states a, where a represents a generalized spin or 
angular momentum, a —] or [. The trapping potential 
is assumed to be harmonic and axially symmetric, 
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Here, pi = ^/xf 



dial and longitudinal coordinate of the ith atom; LOp 



and Zi denote, respectively, the ra- 

and 

uj^ denote, respectively, the angular frequency in the ra- 
dial and longitudinal direction; and iV denotes the total 
number of atoms. We require that the anisotropy pa- 
rameter A ^ 1, where oj^ = Awp, so that the transverse 
motion is "frozen" to zero point oscillations. At zero 
temperature this implies that the Fermi energy associ- 
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FIG. 1: Effective Id coupling constant gid [solid line, Eq. @], 
together with effective Id scattering length aid [dashed line, 
Eq. @] as a function of a^d- 



FIG. 2: Ehom/N (solid line), ^hom (dashed line) and c (inset) 
for a homogeneous two-component Id Fermi g func- 
tion of 7 (horizontal arrows indicate the asymptotic values of 
Ehom/N, Hhom and c, respectively). 



ated with the longitudinal motion of the atoms in the ab- 
sence of interactions, ep = NTiiOz/'^', is much smaller than 
the separation between the levels in the transverse direc- 
tion, tp ^ hcup. This condition is fulfilled if NX <C 1. 
The outlined scenario can be realized experimentally with 
present-day technology using optical traps. 

If the Fermi gas is kinematically in Id, it can be de- 
scribed by an effective Id Hamiltonian with contact in- 
teractions, 



H = NflUn 
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where 



fi2 ^ f)-2 



9idYY5{zi~ zj) (3) 

i=i j=i 



and iV = iV| + iV| . This effective Hamiltonian accounts 
for the interspecies atom-atom interactions, which are 
parameterized by the 3d s-wave scattering length asd, 
through the effective Id coupling constant gid 



gid 



1 



ma^ 1 - Aa^d/ap 



(4) 



but neglects the typically much weaker intraspecies p- 
wave interactions. In Eq. Op = y^h/mujp is the char- 
acteristic oscillator length in the transverse direction and 
A = |C(l/2)|/%/2 ~ 1.0326. Alternatively, gu can be ex- 
pressed through the effective Id scattering length aid, 
gid = -2?i^/(maid), where 



aid 
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Figure ^ shows gid and aid as a function of the 3d s- 
wave scattering length 03^, which can be varied contin- 
uously by application of an external field. The effective 



Id interaction is repulsive, gid > 0, for < 03^ < a^j^ 
(flg^ — 0.9684ap), and attractive, gid < 0, for 03^ > a'^j^ 
and for 03^ < 0. By varying 03^, it is possible to 
go adiabatically from the strongly-interacting repulsive 
regime (gid —>■ +00 or 03^ ^ o,^d)^ through the weakly- 



interacting regime (gid 
attractive regime [gid — 



0) to the strongly-interacting 



-00 or asd 



> 



^3d 



For two fermions with different spin the Hamiltonian 
Eq. ||2J), supports one bound state with binding 



-"id' 



energy ebound 



'h / (mafd) and spatial extent 



aid 

for gid < 0, and no bound state for gid > 0, that is, 
the molecular state becomes exceedingly weakly-bound 
and spatially-delocalized as gid ^ 0~. In the following 
we investigate the properties of a gas with N fermions, 
— 7V|, for both effectively attractive and repulsive Id 
interactions with and without longitudinal confinement. 

To start with, consider the Hamiltonian H^d' ^'i- ®i 
which describes a homogeneous Id two-component Fermi 
gas. The ground state energy Ehom of 77°^ has been cal- 
culated exactly using Bethe's ansatz for both attractive 
|llj| and repulsive [13] interactions, and can be expressed 
in terms of the linear number density uid, nid — N/L, 
where L is the size of the system, 



hom 



Id 



N 



2m 



e(7) 



(6) 



The dimensionless parameter 7 is proportional to the 
coupling constant gid, 7 = mgid/ {h^nid), while its ab- 
solute value is inversely proportional to the Id gas pa- 
rameter nirfjaidj, I7I = 2/nid\aid\- The function 6(7) is 
obtained by solving a set of integral equations llTll . which 
is similar to that derived by Lieb and Liniger for Id 
bosons with repulsive contact interactions. To obtain the 
energy per particle, Eq. 10, we solve these integral equa- 
tions for attractive (7 < 0) jlj and repulsive interactions 

(7 > 0) 111. 

Figure 121 shows the energy per particle, Ehom/N (solid 
line), the chemical potential ^hom, fJ-hom = dEhomI dN 



3 



(dashed line), and the velocity of sound c (inset), which 
is obtained from the compressibility of the system, mc^ = 
niddfihom/dnid, as a function of the interaction strength 
7. In the weak coupling limit, I7I <^ 1, fihom is given by 



TT 



2m 



+ 7 



'Id 



2m 



(7) 



where the first term on the right hand side is the energy 
of an ideal two-component atomic Fermi gas, and the se- 
cond term is the mean-field energy, which accounts for in- 
teractions. The chemical potential increases with increas- 
ing 7, and reaches an asymptotic value for 7 ^ 00 [lo| 
(indicated by a horizontal arrow in Fig.|21), 
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The first term on the right hand side coincides with the 
chemical potential of a one-component ideal Id Fermi gas 
with atoms. Interestingly, for 7^1, the strong atom- 
atom repulsion between atoms with different spin plays 
the role of an effective Pauli principle. 

For attractive interactions and large enough I7I the 
energy per particle is negative (see Fig. I^J, reflecting the 
existence of a molecular Bose gas, which consists of N/2 
diatomic molecules with binding energy Cbound- Each 
molecule is comprised of two atoms with different spin. 
In the limit 7 — s- — 00, the chemical potential becomes 
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The first term is simply tbound/'^, one half of the bind- 
ing energy of the Id molecule, while the second term is 
identical to the chemical potential of a bosonic Tonks- 
Girardeau gas with density nid/2, consisting of N/2 
molecules with mass 2m [l9| . Importantly, the compress- 
ibility remains positive for 7 —00 [a horizontal arrow 
in the inset of Fig. |21 indicates the asymptotic value of c, 
c = Trhnid/ (im)], which implies that two-component Id 
Fermi gases are mechanically stable even in the strongly- 
attractive regime. In contrast, the ground state of Id 
Bose gases with gid < has negative compressibility [2^ 
and is hence mechanically unstable. The formation of 
a quasi- Id molecular gas discussed here has some analo- 
gies to the formation of a 3d molecular BEG |2 1 | and the 
BGS-BEG cross-over for 3d Fermi systems \W^ . 

Using the solutions for the homogeneous two- 
component Id Fermi gas, we now describe the inhomo- 
geneous gas, Eq. within the LDA H H 113. This 
approximation is applicable if the size R of the cloud 
is much larger than the harmonic oscillator length 
in the longitudinal direction, = y^Ti/mUz, implying 
ep ^ fiLdz and ^ 1. The chemical potential of the 
inhomogcncous system can be determined from the local 
equilibrium condition. 




FIG. 3: Energy per particle, E /N (solid lines), and size of the 
cloud, R (dashed lines), for an inhomogeneous two component 
Id Fermi gas as a function of the coupling strength Na^^/a?^ 
for repulsive [gid > 0) and attractive {g\d < 0) interactions. 



and the normalization condition N — J_j^nid{z)dz, 
where z is measured from the center of the trap and 
R — \j2[i' /{muj^^ where [}! = [x for g\d > and 
11' = 11 + \(-hound\/'2 for g\d < 0. The normalization con- 
dition can be reexpressed in terms of the dimensionless 
chemical potential fl and the dimensionless density fiid 
[jl = p! /{y? /2ma\^ and h^d = |al^^|?^ld], 



N- 



'Id 



dx 



(11) 



This expression emphasizes that the coupling strength 



is determined by Na^^/aj. 



\d/'^'z^ ^^id/^'i ^ 1 corresponds 
to the weak coupling and Naf^/aj, ^ 1 to the strong 
coupling regime, irrespective of whether the interactions 
are attractive or repulsive [Bi]. 

Figure 131 shows the energy per particle E/N and the 
size R of the cloud as a function of the coupling strength 
Na\j^/ aj, for positive and negative gid calculated within 
the LDA for an inhomogeneous two-component Id Fermi 
gas. Gompared to the non-interacting gas, for which 
R — \/Naz, R increases for repulsive interactions and 
decreases for attractive interactions. For Na\j^/a1 <C 1, 
R reaches the asymptotic value ^/2Naz for t he str ongly 
repulsive regime, gid +00, and the value N/2az for 
the strongly attractive regime, gid — > —00. The shrinking 
of the cloud for attractive interactions reflects the forma- 
tion of tightly bound molecules. In the limit gid —00, 
the energy per particle approaches e{,ouTid/2 -I- NhuJz/4,, 
indicating the formation of a molecular bosonic Tonks- 
Girardeau gas, consisting of N/2 molecules. 

Using a sum rule approach, the frequency uj of the 
lowest compressional (breathing) mode of harmonically 
trapped Id gases can be calculated from the mean-square 



size of the cloud (z^) [3, 



M = IJ-homlnidiz)] 
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FIG. 4: Square of the lowest breathing mode frequency, uP" , as 
a function of the coupling strength Nal^/a^ for an inhomoge- 
neous two-component Id Fermi gas with repulsive [gid > 0) 
and attractive {gid < 0) interactions determined numerically 
from Eq. (1121 (solid lines). Dashed lines show analytic ex- 
pansions (see text), which are valid in the weak and strong 
coupling regime. 



In the weak and strong coupling regime {Na\^/a1 3> 1 
and <C 1, respectively), {z^) has the same dependence on 
LOz as the ideal Id Fermi gas. Consequently, lo is in these 
limits given by 2^2, irrespective of whether the interac- 
tion is repulsive or attractive. Solid lines in Fig. |31 show 
up , determined numerically from Eq. (|12|l . as a function 
of the interaction strength Nal^/al- A non-trivial be- 
havior of LU^ as a function of Naf^/al is visible. To gain 



further insight, we calculate the first correction 6uj to the 
breathing mode frequency uj [u> = 2uJz{1 + Suj/tOz + •••)] 
analytically for weak repulsive and attractive interac- 
tions, as well as for strong repulsive and attractive inter- 
actions. For the weak coupling regime, we find Slu/uJz — 
±{4:/37r^) / {Nal^/aiyP , where the minus sign applies 
to repulsive interactions and the plus sign to attractive 
interactions. For the strong coupling regime, we find 
Siv/ujz = -[16%/21n(2)/157r2](7Va2^/a2)i/2 for repulsive 
interactions and Slu/lUz ~ (8A/2/157r^)(iVa^^/a^)^/^ for 
attractive interactions. Dashed lines in Fig. 0] show the 
resulting analytic expansions for w^. These expansions 
describe the lowest breathing mode frequency quite well 
over a fairly large range of interaction strengths. Not 
surprisingly, their validity breaks down for Na^j^/a^ ^ 1. 

In conclusion, we have investigated the cross-over from 
weak to strong coupling of quasi- Id harmonically trapped 
two-component Fermi gases with both repulsive and at- 
tractive effective interactions. The frequency of the low- 
est breathing mode, which can provide an experimental 
signature of the cross-over, is calculated. We predict the 
existence of a stable molecular Tonk- Girardeau gas in 
the strongly attractive regime. We believe that our re- 
sults, which are derived using a Id model Hamiltonian 
with contact interactions, hold when more realistic inter- 
actions are considered and when the transverse motion 
is treated explicitly. 
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